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ABSTRACT 

Several GRB afterglow light curves deviate strongly from the power law decay observed in most bursts. We 
show that these variations can be accounted for by including refreshed shocks in the standard fireball model 
previously used to interpret the overall afterglow behavior. As an example we consider GRB 021004 that 
exhibited strong light curve variations and has a reasonably well time-resolved polarimetry. We show that the 
light curves in the 7?-band, X-rays and in the radio can be accounted for by four energy injection episodes in 
addition to the initial event. The polarization variations are shown to be a consequence of the injections. 
Subject headings: gamma rays: bursts — gamma rays: theory — polarization 



L INTRODUCTION 

A long-duration gamma-ray burst (GRB) is now generally 
believ ed to occur following the core collapse of a massiv e 
star (Woosiev 199?, iStanek et alJ l2()()3': 'Hi orth et al.ll2()0l . 
Following the collapse, the released energy pierces a hole 
through the star along its rotation axis, sweeps up and shocks 
the surrounding interstellar medium and produces a long-lived 
afterglow emission. According to the standard fireball model 
(|Meszaros"2002'; 'Zhang & Meszaros '2003'; ■piran"2004), the 
expanding fireball slows down as it sweeps up more and more 
ambient material and decays in brightness. 

The standard fireball model (SFM) predicts that the after- 
glow light curves should be power laws in time, with a break 
due either to one of the characteristic frequencies of syn- 
chrotron radiation moving through the observing band (e.g. 
|Sari et al. 1998) or because of the outflow being collimated 
iRhoads 1999) . The former would only give rise to a modest 
steepening, while in the latter case it can be substantial in ad- 
dition to being achromatic. In most cases where a break has 
been observed, the steepening is substantial and achromatic 
in the optical domain and the latter explanation appears to be 
the appropriate one. 

For the first GRB afterglow light curves observed, the SFM 
provided adequate interpretation of the data. In recent years, 
several well observed afterglows have shown s trong deviation 
from smooth power-la ws, e.g. GRB 011211 ( Jakobsson et alj 
12004) . GRB 0210 04 (iFox et a l. 2003; Holland et al. 2003^ 
and GRB 030329 (iPrice et al.; 2003; Matheson et al. 2003lJ 
iLipkin et al."2004'). Bursts like GRB 000301C, GRB 020813 
and GRB 010222 also deviate from the standard behavior as 
their initial light curve decay rate was slower than can be eas- 
ily accommodated within the model. In these cases, a contin- 
uous energy injection may e xplain the slow decay a s has been 
suggested for GRB 010222 (>Biorn sson et al.l2002h . 

For light curves showing 'bumps' and 'wiggles' during 
the first few days, as e.g. GRB 021004, it has been sug- 
gested that these might arise when the fire ball encounters den- 
sity irregularities in the ambient medium jLazzati et alJ2002t 
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iHevl & Pernall2003t iNakaret alJ l2003h . Higher density re- 
gions would cause brightening episodes, while a rarification 
immediately following the higher density a ppears necessary 
to get the light curve decay 'back on track'. Fox et alJ J2003h 
suggested that the light curve variability of GRB 021004 may 
be due to refreshed shocks. Recently, it has been claimed that 
such models can be rejected as they are unable to explain si- 
multaneously the bumpy light curves and the polarization as 
the latter would not be affected by the freshly injected energy 
(Lazzati_eLal. 2003, hereafter L03). 

In this Letter we show, contrary to the conclusions of 
L03, that several energy injection episodes may in fact ex- 
plain the afterglow light curve re-brightenings of bursts like 
GRB 021004, each injection contributing to the light curve 
as it catche s up with the pre viously expanding shock front 
(see also Granot et alJ 120031 for a similar conclusion on 
GRB 030329). We consider injections in the SFM and ap- 
ply our model to GRB 021004. We show that not only are 
the light curves readily accounted for, but also that its polar- 
ization properties follow directly from the model. Detailed 
account of this work as well as application to other bursts will 
be presented elsewhere (Johannesson et al. in preparation). 

2. ENERGY INJECTION 

For a gen e ral rev i ew of the SFM w e refer the reader 
to iMeszarosI IZhang & Meszaros! J2003h and iPirad 

i2004l) . The basic idea behind the model is the self-similar 
relativistic shock solution of Blandford & McKee ( 1976\ 
where it is assumed that energy is released instantaneously 
at the onset of the burst. It is, however, both natural and 
expected that the energy release may be either continuous 
or episodic (Rees & Meszaros 1998; Kumar & Piran 200^ 
Sai-i & Meszaros 2000). 

We extend the SFM as suggested by iRees & MeszarosI 
Jl998t by applyi ng energy and momentum conservation, as 
in Rhoads ( 1999), to both discrete and continuous energy in- 
jections although only discrete injections will be considered 
here. We assume that several discrete shells are injected si- 
multaneously with different Lorentz factors. The shell with 
the highest Lorentz factor, Fq, drives the initial evolution of 
the afterglow. Once it has decelerated to a value lower than 
the Lorentz factor of the next shell, they will collide with a de- 
lay corresponding to the time it takes the second shell to reach 
the first. The collision is assumed to be instantaneous and its 
dynamics is neglected as is any anisotropy that may results 
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from the shell interactions. In afterglow interpretations, the 
delay is fixed by the start of a brightening episode. The en- 
ergy of the shell is determined by the observed increase in flux 
level. Similar results can be ac hieved with a continuous en- 
ergy injection as discussed by Panaites cu. Meszaros & Reesi 
Iil998 ). but requires a very steep energy injection profile and 
in some cases even a varying electron energy index, p. 

In calculating light curves and spectra, we assume that the 
radiation is of synchrotron origin, and we consider the lo- 
cal synchrotron spectra at each point in the outflow to con- 
sist of power law segments, smoothly joined at the charac- 
teristic frequencies (Granot & Sari 2002). We integrate over 
a thin shell at the equal arrival time surface. Each shell el- 
ement is assumed to be locally homogeneous, its thickness 
being determined by the jump conditions across the shock 
(iBlandford & McKee 1976) and the conservation of particles. 
We can therefore consider general density profiles, for exam- 
ple a constant density environment, a wind or various density 
irregularities. 

We calculate the instantaneous fireball polarization as in 
iRossi et al. ( 2004). We assume a random magnetic field com- 
pressed by the blast-wave, thereby introducin g some align - 
ment perpendicular to the compressed direction (^ Laingl 1 9801) . 
We evaluate at a given observer time the contribution of each 
surface element of the equal arrival time surface to the Stokes 
parameters dU = P{9)dLcos{2(l)), and dQ = P(9)dLsm(2(j)), 
with the angular dependence of the polarization given by 
P(e') = Posin^6'7(l+cos2e'') lLaini lf!980). Here, 9' is the 
angle from the line of sight in the comoving frame and dL 
is the local luminosity. In the optical range for most of 
the fireball evolution, the maximum degree of polarization 
is taken to be Po = (p+ l)/ (p + 7/3) w 70% for p = 2.2 (see 
iRvbicki & LightmanI 19791 for details). For a measurable po- 
larization to occur, the line of sig ht to the observer has to 
be off the outflow symmetry axis ( Ghisellini & Lazzati 1991 
|£ari 1999). As in Ghis elUni & Lazz ati (1999), we find that 
for an instantaneous energy release, the polarization light 
curve has two extrema, bracketing the time when 6* « l/F. 
At ap proximately that tim e the polar izatio n angle rotates by 
90° MGhisel Uni & Lazzati 1999; Sari 1999: Gran ot & Konigll 
l2003h . It is important to note that the maximum degree of po- 
larization and time of change of the polarization angle depend 
strongly on the viewing angle. The polarization evolution de- 
pends mostly on the evolution of F, but is also af f ected by the 
evolution of the jet opening angle. R ossi et alJ (|2004) have 
shown that the larger the lateral velocity, the lower the polar- 
ization. We assume the sideways expansion is given by the 
comoving sound speed, c, = c/^/3. Recent work has shown 
that the lateral expansion may in fact be slower ( Granot et a^ 
[2001 ), but we adopt this assumption to highlight the effects of 
the energy injections. 

The number of model parameters can be quite large if there 
are many brightenings. The global parameters: initial en- 
ergy, £0, Lorentz factor, Fq, half-opening angle, 6*0, ambi- 
ent density, «o (or density profile), relative energy density 
in electrons, eg, and magnetic field, eg, and electron index, 
p, are determined, as in the SFM, from the initial afterglow 
evolution and the total flux level (using the burst redshift). 
All episode parameters are fixed by the observed time of 
brightenings and the increase in flux levels. As discussed by 
|Zhang & Meszaros (2002) (hereafter ZM02), only the energy 
ratios and the relative velocities are relevant as long as the rel- 
ativistic phase lasts, and therefore absolute shell energies and 



Lorentz factors are not needed. If polarimetry is available this 
may be the most reliable way to determine the so called jet 
break time as the polarization angle is predicted to change by 
90° at that time. Determining the break time from a broken 
power-law fit to the optical light curve can lead to an erro- 
neous result if the light curve is not smooth, as e.g. in the case 
of GRB 021004. In addition, for the interpretation to be con- 
sistent, the polarization levels predicted by the model should 
agree with those observed, once all parameters have been de- 
termined. Our code has been extensively tested and compared 
to other work. We have also compared the output with ana- 
lytical results, and find good agreement in parameter ranges 
where the assumptions used in analytical work are valid. 

The effects of a single injection episode can be described as 
follows: When a relativistic shell catches up with the slower 
shock front propagating into the ambient medium, it increases 
the F of the forward shock (ZM02), but the subsequent evo- 
lution of F continues with the same decay rate as before. Re- 
verse shocks may be expected, but these are expected to be 
weak for mild energy injections (ZM02). The flux from the 
fireball increases, but as the energy addition is instantaneous 
the light curve decay will from then on also continue with the 
same rate as prior to the shell injection. The rise in flux will 
not be sharply defined in time as the flux is obtained by in- 
tegration over the equal arrival time surface that smooths out 
the transition. The injected shell will thus result in a smooth 
bump in the light curve. The net result is that the light curve 
is 'shifted' upwards at the time of the injection, but retains it 
original decay rate behavior from then on, the rate being most 
strongly determined by p, and the density structure of the 
ambient medium (e.g. Panaitescu, Meszaros & Rees 1998). 
Each energy injection episode delays the decaying evolution 
of F, and causes it to increase temporarily. In addition, the 
emitting region of the relativistic outflow centered on the line 
of sight, temporarily brightens up and outshines the bright 
emitting ring-like region around the line of sight (see Waxma^ 
1997; Panai tescu & Mesza ros 1998; Granot et al. 1999, for a 
detailed discussion). As a result the flux will increase slowing 
the light curve decay ( Kumar & Piran 200(|). It is important to 
realize that contrary to earlier statements (L03), polarization 
will also be affected by the energy injection, as increasing F 
increases the flux, causes increased aberration, decreases the 
emitting surface area and thus decreases the degree of polar- 
ization compared to an unperturbed evolution. We will show 
examples of these effects in the next section. Additional in- 
jections can be viewed as superpositions of repeated single 
episodes with similar effects on the light curve. 

3. GRB 021004 

As an example, we consider GRB 021004. It showed strong 
light curve variations with a best fit light curve break time 
of 4.74 days (iHolland et al.ll2003l) . This break time was ob- 
tained by fitting a broken power law to the data and resulted 
in a rather large formal error. The afterglow also exhibited 
variable polarization levels and a 90° change in polarization 
angle at approximately 1.0 day (Rol et al. 2003). These two 
time estimates, if interpreted within the SFM, should be sim- 
ilar but are in this case inconsistent with each other. 

We find that the light curves can be explained by 4 injection 
episodes, superimposed on the initial GRB event. We sum- 
marize the model parameters in Tabled and show the model 
light curves in Fig. We set p = 2.2, as theoretical studies 
suggest this to be a universal value (see Piran 2004, for a re- 
view and references). Using the burst redshift, z = 2.335 (e.g. 
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TABLE 1 

Model parameters for GRB 021004 



Parameter 


Value 


£o 


1.0 


To 


800 


P 


2.2 


no 


26.0 


El 


3.5 


Ez 


5.6 


£3 


13.0 


£4 


7.0 



Initial half-opening angle is = 1.4°, the line of sight angle is 9y = 0.95^?o• 
tg = 0.21, and eg = 2 ■ 10"^. «o is in units of cm"^, £0 is in units of 10^" ergs, 
and the four energy injection values, £",, are relative to Eq. The four injection 
times in the observer frame are approximately 1 h, 16h, 42h and 105 h. 
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Light curves of GRB 021004: 


i?-band (dotted), 8.^ 



(solid), 22.5 GHz (dashed), and X-ray (multiplied by a factor of 100; dash- 
dotted). Variability in the optical light curve is obtained by imposing 4 en- 
ergy injection episodes at times indicated by the arrows. Parameters were 
adjusted so that the model Hght curves would also go through the X-ray 
and radio light curves as well as the polarization data. A host galaxy 
of R-band magnit ude 24 was a dded to the model optical light c urve. R- 
band d ata is from Uemura et al. i2003); Pandev et al. i2003); Holland et alj 
120031). X-ray data is from Sako & Hamson 1 2002a b). Radio data is from 
IBerger. Frail & Kulkarni 1 2002); Frail & Berger 1 2002 ) and the GRB Large 
Program at the VLA (AK509)' . Radio data hints at a 5th injection at 1 1-12 
days. 



lM0ller et alJ200l . Eq, Tq and «(), are determined by demand- 
ing that the model flux matches that observed in the same way 
as in the standard model. The bumps in the /?-band light curve 
are then used to set the times and relative energies of the in- 
jection episodes. Other parameters, such as e^, and e^, are then 
adjusted until the model agrees with the radio and X-ray data. 
Finally, the change in polarization angle to fixes the opening 
angle of the jet, 1.4°. Recall that this depends on the rate of 
lateral expansion, here assumed to be constant. 

All four injections are mild (ZM02), the relative Lorentz 
factor in all cases about 2. The relative energy of the first 
injection is 3.5, in subsequent episodes it is 1.24, 1.27 and 
about 0.3 of the total in the fourth episode. The peak flux from 
the reverse shocks is maximum for the first event, a factor 
of 30 higher than from the forward shock with the frequency 
at maximum reverse flux of about 2.5 GHz. No radio data 
is available at the time of the first injection. All subsequent 
reverse shocks have lower maximum fluxes and frequencies in 
the sub-GHz range. We therefore neglected the reverse shock 
contributions to the flux. 
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Fig. 2. — A segment of the R-band light curve of GRB 021004 (top panel). 
Effects of each of the energy injection events is shown. Dotted curve shows 
the theoretical afterglow model without additional energy injection. Dashed 
curve shows the effect of one injection event, dot-dashed two events, dash- 
triple-dotted three events and the solid curve shows the effect of all four in- 
jection events. Arrows indicate the injection times. Lower panel shows the 
corresponding polarization light curves. Note how each injection episode 
causes the polarization to deviate from the curve of unperturbed evolution. 
Polarimetry data is from .Rol et al. 1 2Q03J and L03. 



As seen in Fig. the model is able to account for the 
light curves in all observed wavelength regions at all times, 
except for the first half hour or s o when an initial reverse 
shock may be dominating the flux ( Kobavashi & Zhangl2003t 
Panaitescu & Kumar 2004), and the radio data after about 20 
days. A fifth injection may be able to account for the late 
radio data followed by a transition to a non-relativistic expan- 
sion regime at 50-60 days. 

The intrinsic optical and X-ray spectral slopes are fixed by 
our choice of p = 2.2. At 1.4 days these are /3opt = l)/2 = 
0.6 and /3x = /'/2 = 1.1. The former is within 2a of the 
intrinsic (e xtinction corrected) spectral index, estimat ed in 
iHolland et a l. (2003), of Pun = 0.39 ±0.12. Matheso n et all 
ll2003ah find a steeper slope of the optical spectra but with a 
clear curvature, (3 = 0.98 ± 0.03. [Holland et al. (2003) also 
obtained an excellent fit to the 2-10 keV Chandra spectral 
index of 1.03 ±0.06, agreeing with the model prediction. 
In the 0.4-10 keV range and including absorption, they find 
Px = 0.94 ±0.03. These results are consistent with those of 
Butler et al. (2003) and Sako & Harrison (2002a). 

In Fig.|2]we plot the /?-band and polarization light curves 
of GRB 021004. We show the effect of each injection episode 
separately as well as the combine d result. W ith this inter- 
pretation, the estimate bv Holland et alJ J2003i) of the late jet 
break time now has a simple explanation. The repeated en- 
ergy injections slow the early light curve decay and delay the 
steepening until after the last injection. A broken power law 
thus underestimates the pre-break decay slope and overesti- 
mates the jet break time. The change in polarization angle 
should give a more reliable estimate of the break time. 

As mentioned above, with only one energy release event, 
we rec over the polarization behavior of Ghisellini & Lazzay 
( 1999) with differences due to our more detailed fireball 
model (dotted curves). Adding one injection episode clearly 
shows the effects on the optical and polarization light curves 
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(dashed in fig. |2ji. At the time of injection the polarization 
level drops sharply and reaches a minimum just after max- 
imum brightness in the light curve. In this particular burst, 
the first injection occurs before the jet break, and because the 
injection modifies the temporal evolution of F, it delays the 
jet break time compared to a single event, and hence also the 
90° change in polarization angle. The break time, as defined 
by the change in polarization angle, is approximately at 0.6 
days. This is just before the 2nd injection and therefore goes 
unnoticed in the light curves until after the last injection. 

Including all four injections results in the polarization light 
curve shown in the lower panel in fig.|2l(solid). There, the ef- 
fect of each injection episode is very clear Depending on the 
strength of the injections, the polarization angle may change 
again, although at a very low polarization level (see local min- 
imum at approximately 4 days). The correlation between the 
observed flux and polarization variations is in this interpreta- 
tion seen to be directly rooted in the dynamics of the outflow. 

4. DISCUSSION 

We have also applied our code to a SFM modified with den- 
sity variations in a homogeneous medium, but without energy 
injections. We modeled the variations with Gaussian profiles 
as in L03. The calculated flux is very sensitive to the num- 
ber of radiating electrons and care must be taken in counting 
them. Only when we specifically introduced the shock profile 
of Blandford & McKee ( 1976), did we manage to get suffi- 
cient brightenings. A proper treatment of the effects of den- 
sity variations, valid at all observer times and radial ranges of 
interest here, requires a numerical solution of the dynamical 
equations and a self consistent integration over the emitting 

^Jakar&Orei] (12004 1) have interpreted the GRB 021004 op- 
tical data using the "patchy shell" model where the angular 
distribution of energy in the shock is inhomogeneous due to 
random fluctuations. Its application requires the angular en- 
ergy distribution to be specified, but with suitable parameters 
the model is able to account for the general trends in the data. 

A contribution to the polarization may originate in 



the Galaxy or within the host. L03 conclude that for 
GRB 021004, the host contribution can be accounted for by 
using spectropolarimetric data, while the Galaxy may domi- 
nate at low polarization levels. However, the time variation of 
the polarization can only originate in the source. 

It is natural to assume that the central source releases energy 
in several discrete events, essentially simultaneously. The to- 
tal energy injected into the collimated outflow inferred from 
our model is about 3 • 10^' erg. It is of the same o rder as es- 
timated in other bursts after beaming correction (Frai l et all 
2001). Applying b eaming correctio n to the isotropic 7-ray 
energy estimated in iFox et'aP (l2003h . we find that 1.5 • 10"^' 
ergs were emitted in 7-rays. We remind the reader that the 
narrow opening angle inferred from our modeling is a conse- 
quence of assuming a constant lateral expansion velocity. In 
this interpretation, most of the electromagnetic energy output 
in GRB 021004 was emitted at longer wavelengths. 

In addition to GRB 021004, we have considered other burst 
with highly variable light curves, including GRB 970508 and 
GRB 030329. Most cases considered so far can be accounted 
for by discrete energy injections, with GRB 030329 a possible 
exception. This work will be discussed separately. 

The self similar solution is based on the assumption of a 
uniform or a wind structured ambient medium. It may there- 
fore be able to account for the overall behavior of the expand- 
ing shock, but one cannot expect it to be able to follow den- 
sity variations in detail. For that a numerical solution of the 
dynamical equations is required. We conclude that refreshed 
shocks can in fact account for the observed variability both 
in the light curves and in the polarization properties of bursts. 
The strongest argument in favor of our interpretation is the 
fact that we are able to account for both broad band behavior 
as well as the polarimetry within a single model. 
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